Landau levels and oscillator strength in a biased bilayer of graphene 
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We obtain analytical expressions for the eigenstates and the Landau level spectrum of biased 
graphene bilayers in a magnetic field. The calculations are performed in the context of a four- 
band continuum model and generalize previous approximate results. Solutions are presented for the 
spectrum as a function of interlayer coupling, the potential difference between the layers and the 
magnetic field. The explicit expressions allow us to calculate the oscillator strength and the selection 
rules for electric dipole transitions between the Landau states. Some transitions are significantly 
shifted in energy relative to those in an unbiased bialyer and exhibit a very different magnetic field 
dependence. 
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I. INTRODUCTION 

The recent realization of stable single layer and bilayer 
carbon crystals (graphene) has aroused considerable in- 
terest in the study of their electronic propertiesii^ii^ii^i^. 
These materials have been shown to display interest- 
ing new phenomena that may allow the development of 
carbon-based nanoelectronic devices. The behavior of 
charge carriers in wide single-layer graphene sheets has 
been described as that of ultrarelativistic massless parti- 
cles, with a "light speed" equal to the Fermi velocity of 
the crystal and a gapless linear dispersion close to the K 
and K' points. One consequence of that is that single- 
layer graphene displays an unusual quantum Hall effect, 
in which the quantum Hall plateaus are found in half- 
integer multiples oiAe^/h. Moreover, the massless char- 
acter of the charge carriers in single-layer graphene gives 
a dependence to the Landau levels (LL), which has 
been recently confirmed by infrared transmission and cy- 
clotron resonance experiments in thin graphite samples^ 
and on single layers of graphene^. On the other hand, for 
an unbiased graphene bilayer the spectrum at the vicin- 
ity of the K points displays four parabolic bands but is 
still gapless^iiS. The absence of a gap, together with the 
chiral nature of the electronic states, in both single-layer 
and bilayer graphene, is at the root of phenomena such 
as the Klein tunnelin g^^i^^i^^ which has important con- 
sequences for the design of graphene-based devicesii.. 

Recent theoretical and experimental results have 
shown that a gap can be induced in bilayer graphene 
by changing the density of charge carriers in the layers 
through the application of an external field or by chem- 
ical doping, which creates a potential difference between 
the layersi^ii^. The presence of such a bias transforms 
the graphene bilayer into a semiconductor with a tun- 
able gap. This implies that such biased bilayers can 
be used to confine charge carriers through the use of, 
e.g., nanostructured gates or spatially varying doping^^. 
Tight-binding calculations have shown that the existence 
of such a gap can have a significant effect on the LL spec- 
trum of biased graphene bilayers22i^ 



In this work we determine the spectrum of a biased 
graphene bilayer in the presence of an external magnetic 
field in the continuum description and obtain, by solv- 
ing the four-band Hamiltonian, the exact eigenstates as 
well as the oscillator strengths and the selection rules for 
the electric dipole transitions between LL. Previous nu- 
merical results for the LL in biased and unbiased bilayer 
graphene2ii2^, were obtained within a model based on a 
reduced two-band Hamiltonian, which is accurate for en- 
ergies much lower than the interlayer coupling parameter. 
Furthermore, these previous studies used a description in 
which the Hamiltonian can be described approximately 
in terms of ladder operators in the basis of Landau func- 
tions, which is only possible in the Landau gauge. Here, 
we generalize these results to much higher energy by us- 
ing the more accurate four-band model, and obtain exact 
expressions in both the Landau gauge and the symmetric 
gauge. The calculations allow us to obtain a significant 
correction to the spectrum obtained by previous approx- 
imate methods. The model is presented in Sec. H for 
zero magnetic field B. A finite field B is considered in 
Sec. HI, in the Landau and symmetric gauges, and the 
results are given in Sec. IV. The oscillator strength is 
evaluated in Sec. V and a summary follows in Sec. VI. 



II. ZERO MAGNETIC FIELD 

The crystal structure of an undoped bilayer of 
graphene is that of two honeycomb sheets of covalent- 
bond carbon atoms coupled by weak Van der Waals 
forces. To each carbon atom corresponds a valence elec- 
tron and the structure can be described in terms of four 
sublattices, labelled A, B (upper layer) and A', B' (lower 
layer). In this work we focus on the Bernal stacking, in 
which half the atoms of the upper layer (i.e. the sites of 
the A sublattice) are on top of half the sites of the lower 
layer (the sites of the B' sublattice). This type of stack- 
ing is found in graphite and has been reported as well for 
samples of bilayers of graphene. Alternatively, one can 
consider the hexagonal stacking, in which all sites of the 
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top layer are on top of all sites of the lower layer. That 
configuration introduces a different coupling scheme be- 
tween the layers, and can lead to qualitatively distinct 
results. However, the band structure data from the epi- 
taxial bilayer samples of ref. 15, as well as the results 
from the samples produced by micro-mechanical cleav- 
age of graphite in ref. 21 are consistent with the Bernal 
stacking. The coupling between the layers is described 
by an interaction term between the A and B' sublattice 
sites. Considering only nearest-neighbor hopping, the 
Hamiltonian of the system in the vicinity of the K point 
is given, in the continuum approximation, by^l 
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H = Ho + (A[//2)t, 



with 
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the operator Tz assigns a positive (negative) sign to the 
upper (lower) layer labels and is defined as 



I 

-I 
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with I denoting the 2x2 identity matrix. We assume 
a constant interlayer coupling term t ^ 400 meV; tt — 
vpiPx+iPy), P = (Px^Py) is the 2D momentum operator, 
vp^lx 10^ m/s, Uo = {Ui + U2)/2, AU = Ui- U2 and 
Ui, U2 are the potentials at the two layers, which reflect 
the influence of doping on one of them and/or the inter- 
action with an external electric field, i.e., the gating. In 
this work we treat the potentials at each layer and the gap 
as adjustable parameters, which may be obtained exper- 
imentally, e.g. from angle-resolved photoemission spec- 
troscopy measurementesl^. The eigenstates of Eq. (1) 
are four-component spinors ^' — [tpA , V's , ipB' , i'A']'^ , 
where tPa,b (V'A'.b') are the envelope functions associ- 
ated with the probability amplitudes at the respective 
sublattice sites of the upper (lower) graphene sheet. The 
superscript T denotes the transpose of the [...] vector. 

If the magnetic field is absent and the potentials Ui and 
U2 are constant, the single-particle spectrum consists of 
four bands with eigenvalues^ 



E+{k) = Uo + {l/2)[{t±T)^ + n] 
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FIG. 1: Low-energy electronic dispersion of graphene bilayers 
for B = for different values of the gap, namely AU = 
(black solid line), AU = 25 meV (red dashed line), AU = 100 
meV (blue dotted) and AU = 200 meV (green dot-dashed 
line). In all cases Uo ~ 0. 



line), AU = 100 meV (blue dotted) and AU = 200 
meV (green dot-dashed line). For large fc- values the lin- 
ear E — k behavior is recovered. The hole spectrum 
in this case is obtained by taking E{k) —E(k). In 
the absence of bias, the spectrum shows approximately 
parabolic bands, which become increasingly deformed as 
AU increases, with the appearance of energy minima 
at non-zero values of fc, and can be approximated by a 
fourth-order polynomial function^' As we show below, 
this non-parabolic dispersion has a significant influence 
on the LL results. 



III. FINITE MAGNETIC FIELD 

We consider an undoped bilayer of graphene in the 
presence of an external perpendicular magnetic field B. 
It is instructive to consider separately the two gauges for 
the vector potential A given below. The momentum op- 
erator p, appearing in the term TLa in Eq. (1), is shifted 
by — eA. 



A. Landau gauge 



E±{k) = Uo-{i/2)[{t±rf + n] 
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where T = [t'^ + 43% + A{s%/f)AU^Y' " , sp = hvpk, and 
= [1 - 4s|,/f2] At/2. Note that for fc = the spectrum 
shows a gap at A; = of size E^{0) — EZ{0) = \AU\ and 
the system becomes a narrow-gap semiconductor. Figure 
1 shows the low energy electronic spectrum of graphene 
bilayers for Uq = and three different values of gap: 
AU = (black sohd hue), AU = 25 meV (red dashed 



In this gauge we have A = {0,Bx,0) with [Ho,Py] — 
0. Consequently we can write the solutions for the 
spinor components in the form ipA{x,y) = ^A(a;)e*'^''^, 
tpBix,y) = i(t)B{x)e'^^y , tpA'ix,y) = i(j)A'{x)e'^yy and 
ipB'{x, y) = (/)B'(a:)e*'^»^. The x dependence of the spinor 
components is then described by 



[-^ + ky + f3x^c^A = {E-Ui) 



(6) 
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f(^^ - ky - (3x](j)B ^ -{E - Ui)c^A + t(t)B', (7) 
d 



(^— + ky + f3x'^(l)A' = -{E-U2)(f>B'+t<l)A, (8) 



(9) 



where [3 = eB/h. For a wide, uniform system these equa- 
tions can be decoupled to obtain, for ipA 

{-^ + ^' + ^'- ^'Y^A = '^')*" + [1 - 2Hn0A 

(10) 

with ^ = ^Bky + x/ Ib] the characteristic length scale 
of the system is the magnetic length — [/i/eS]^/^, 
a = e — Mo is the energy shifted by the average potential 
between the layers uq = (iti-|-M2)/2, and 5 = Am/2, Au = 
Ml — M2- The energy, the potentials, and the interlayer 
coupling strength are written in dimensionless units as 
e = EiB/fivp, ui,2 = Ui,2(-B /fi-VF and t' — tiB/hvp, 
respectively. Next, we we consider the equation 



d^/de+e]c^A = i±{a)<j^A, 



(11) 



where the potential-dependent eigenvalues are 

7±(a) = a'^5^± [{a" - 5^)t'^ + (1 - 25af]"\ (12) 

With the substitution (j)A — /a exp(— ^^/2), and after 
some straightforward algebra, the LL spectrum is given 
by the solutions a„ of the fourth-order algebraic equation 

[(a + (5)2-2(n+l)][(a-(5)2-2n]-(a2-,52)i'2 ^ 0, (13) 

where n is an integer. The solutions of Eq. (11) can be 
obtained in terms of Hermite polynomials. 
Substituting (j)A in Eqs. (6)-(9) we obtain 

4>A = dnC„ HniO e~^^/^, 

V2n 



' {an - S) 



c„_ii/„_i(C)e-« 



(an + 0) 

(/-B' =/„dnC„i7„(0e-«'/2, (14) 
where -ffn(C) is the Hermite polynomial of order n and 
c„ = l/(n!2"V^)'/', 
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Notice that the eigenstates are determined by the quan- 
tum numbers (ky , n) while the spectrum is independent 
of k„. 



B. Symmetric gauge 

In problems with cylindrical symmetry it is advan- 
tageous to use the symmetric gauge A = (0, Sp/2,0). 
In this case the corresponding solutions for the four- 
component spinors VE" are written as^^ii^ 



/ 0^(^)6™" \ 

i(/)B(p)e*(™-i)'' 
0i3'(p)e^™^ 
\^i0A'(p)e'("+i)'*y 
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where m is the angular momentum label. The radial 
dependence of the spinor components is described, in di- 
mensionless units, by 
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(pA ^ (a~ S)(j)B, 
- ? 0s = -(a - S)(t)A + t'(j)B', 
+ £, <Pa' = -(a + S)(pB' +t'<l>A, 
'B' = (a + S)(l)A', (17) 



where now ^ — \/2pj (21b\ and the other terms are de- 
fined as before. 

Following the same procedure as in the previous case, 
one can obtain an equation that is identical to that of an 
harmonic oscillator 
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0A = 27±(a)0A, (18) 



where 



l±(a) = (n' + l/2), 



(19) 



with n' = n + (|to,| -|- ni)/2, n — 0,1,2... and 7±(a) is 
given by Eq. (12). The LL spectrum is determined by 

[(a+5f-2(n' + l)][(a~5f~2n']-(a''~5'')t'^^Q, (20) 

which is identical to Eq. (13) except that here the Lan- 
dau index n' is used instead of n. 

It is evident, from Eq. (18) that the spinor components 
(j)A and (pB' can be obtained in terms of generalized La- 
guerre polynomials. For to > we obtain: 
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For m < the result is 



(21) 



r<m+l ^|m|+l -t- ^ |m| + l /^2\ 



\/2(n' + 1) „,^i _^ |m|-1.^2N 

'^B'=C;,\e''"'e^^4™'(^'), (22) 
where 

C;r„ = rfn/ [«!/("^ + |m|)!]^/' , C^,„ = /„C^„; (23) 
the other constants are as before. 

C. Special cases 

i) t' ~ 5 = 0. This is the case of a single graphene 
layer or two uncoupled, unbiased layers. The solutions 
are then 

an' = ±v/2(n' + l), an' = ±V2^, (24) 

which are the well known expressions for the LL in a 
single-layer of graphene. 

ii) 6 = 0. That is the case of an unbiased bilayer. The 
explicit solutions are 
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This expression can be simplified by expanding the inter- 
nal square root. By taking the negative sign, one obtains 



an' =±-y/n'{n' + 1)/ 



l + -(2n' + l) 
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(26) 



which is valid for n'/t''^ << 1. By assuming 1/[1 
(2/i'2)(2n' + l)]i/2 ~ 1, E2. (29) is simplified to 
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n' = ±-V"'(n' + l), 



(27) 



which is the expression previously obtained in the con- 
text of the two-band continuum model^-. By taking the 
positive sign in front of the internal square root in Eq. 
(25), one finds the higher-energy LL that arise from the 
upper band at E — t. For n' = that gives 
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FIG. 2: Landau levels (n = — 5) in an unbiased graphene 
bilayer as a function of the magnetic field. The solid, red lines 
result from a numerical solution of Eq. (25) and the dashed, 
blue ones are the approximate results from Eq. (27). 



IV. NUMERICAL RESULTS FOR THE 
ENERGY SPECTRUM 

Figure 2 shows the first 6 positive energy LL as a func- 
tion of the magnetic field B for an unbiased grahene bi- 
layer. The figure shows that as B increases, Eq. (27) 
consistently overestimates the values of the energy levels 
(dashed lines) as compared to the exact results obtained 
from Eq. (25) (solid) (which coincide with those from the 
approximate expression Eq. (26)). The discrepancy in- 
creases with B and, for the parameters usually observed 
experimentally {t = 0.4 eV, vp ~ 1-0 x 10^ m/s), can at- 
tain 5 meV for B = 15T and n = 1. This result was con- 
firmed by a numerical calculation based on a discretiza- 
tion of Eq. (17). Recent experimental studies have mea- 
sured the LL spectrum in bilayers of graphene for fields 
up to 12 T— . At these fields, the discrepancy between 
the four-band and two-band results are w 10% for the 
n = I LL. This is below the energy resolution of the 
current experimental techniques. However, for stronger 
fields the difference should become large enough to be- 
come detectable. In addition, the presence of the LL 
associated with the higher energy band at E fa 0.4 eV 
may be detectable by the scanning tunneling technique 
discussed in ref. 23. 

Equation (20) shows that the charge-conjugation sym- 
metry between electrons and holes is preserved for the 
Landau states in a single layer {t — 0, 6 = 0) as well as 
in the unbiased double layer (5 = 0), since the equation 
is unchanged under the transformation a — > —a. On the 
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other hand, that ceases to be the case for a biased bilayer 
of graphene. Figure 3 shows the first 6 fow-energy LL, 
for a biased bilayer with C/i = — C/2 = 50 meV. The pos- 
itive (negative) energy levels are shown in Fig. 3a (3b), 
and the results show a non-monotonic dependence on the 
magnetic field. To each index n there correspond four 
solutions, two negative and two positive, with the low- 
energy branches being equal to \6\ for B = 0, whereas 
the high-energy solutions arise from the bands at \a\ — t' 
(not shown). The only exception is the n = state, for 
which, in the case of the K valley, there is no positive 
solution at the vicinity of S (the opposite is true at the 
vicinity of the K' valley). The figure shows the LL cor- 
responding to n = (black solid line) , n — 1 (red dashed 
lines), 2 (blue dotted), 3 (green long-dashed), 4 (orange 
dot-dashed) and n = 5 (magenta dot-dot-dashed lines). 
One striking feature in the spectrum is the appearance 
of crossings and energy minima in some magnetic field 
regions shown in the insets. These features arise as the 
potential difference (i.e. the gap parameter) between the 
layers is increased, with the energy minima and the cross- 
ings initially appearing at very small values of B. As 
AU increases, the minima in each branch ocurr at larger 
magnetic fields and lower energies. These results can be 
explained by taking into account the band structure of 
the biased graphene bilayer. The introduction of a mag- 
netic field in a 2D system constrains the electronic states 
to closed cyclotron orbits, which in turn causes an energy 
quantization in which the energy eigenvalues correspond 
to momentum states that are integer multiples of h/is- 
For a system with a parabolic band this gives rise to a 
linear dependence of the LL on the magnetic field. For 
a single graphene layer, on the other hand, the linear 
momentum dependence of the energy introduces a ^/B 
dependence of the LL spectrum. In a biased graphene 
bilayer the electron dispersion has a more complicated, 
"mexican hat" shape (see Fig. 1) that becomes more 
pronounced as AU becomes larger. For B — this dis- 
persion then allows the existence of pairs of degenerate 
states with different values of momenta which, for finite 
B, can become degenerate LL with different Landau in- 
dices. By setting n' = in Eq. (20) we obtain, for a 
values close to —S, 
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FIG. 3: a) Positive-energy Landau levels in a biased graphene 
bilayer as a function of the magnetic field B. (b) As in (a) 
for negative-energy levels. In both cases Ui — —U2 ~ 50 meV 
and n = (solid line), n = 1 (dashed lines), 2 (dotted), 3 
(long dashed), 4 (dash dotted) and n = 5 (dot-dot-dashed 
lines). The rectangular areas are enlarged in the insets. 
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FIG. 4: Landau levels in an unbiased graphene bilayer, for 
i? = 5 T, as a function of the interlayer coupling parameter t. 
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The introduction of bias, together with the presence of 
a magnetic field, breaks the inversion symmetry of the 
bilayer, which in turn breakes the charge-conjugation 
symmetry. A similar effect is obtained for the case of 
states in a position-dependent confining potentialii. This 
electron-hole asymmetry is particularly evident in the 
lowest LL levels, which are not distributed symmetrically 
around E — Uq, i.e. En ^ —E^n, for Uq = 0, except for 
AU « 0. 

Previous studies of LLs in bilayer graphene that were 
based on a continuum model have used a two-band ap- 
proximation, which is expected to be accurate for large 



coupling interlayer strengths t. The four-band model 
used here allows one to investigate the whole range of 
interlayer interaction strengths. Figure 4 shows the first 
10 LL for an unbiased graphene bilayer as a function of 
the coupling parameter t. The t = limit corresponds 
to the case of two uncoupled single-layers, and the spec- 
trum is given by Eq. (24). As the coupling increases, 
the two-fold degeneracy of the levels is lifted, with the 
higher energy levels being shifted towards E — t, and 
several crossings occur. The figure also shows that the 
energy difference between the lowest LL decreases for in- 
creasing t, with the n = LL remaining unaffected by 
the coupling, in the absence of bias. 

The dependence of the spectrum on the energy gap can 
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FIG. 5: Landau levels in a graphene bilayer as a function of 
the difference AU in potential between the layers for B = 5 
T: n = (black solid line), n — 1 (red dashed lines), 2 (blue 
dotted lines), 3 (green long-dashed lines), 4 (orange dash- 
dotted lines), and n = 5 (magenta dot-dot-dashed lines). 



be seen in Fig. 5, which shows the LL corresponding to 
n = (black solid line), n = 1 (red dashed lines), 2 (blue 
dotted), 3 (green long dashed), 4 (orange dash dotted) 
and n = 5 (magenta dot-dot-dashed lines). One striking 
feature of this result is the reduction of the energy differ- 
ence between some of the LL as the gap increases, with 
the eventual appearance of degeneracies. That can be 
explained by considering the fact that as AU increases, 
the band structure of the bilayer evolves from an approxi- 
mately parabolic dispersion to the approximately fourth- 
order polynomial dependence on the momentum. From 
Eq. (20) one can easily show that degeneracies will occur 
between LL with indices ni and ri2 for energies and fields 
that satisfy the relation 



ni +n2 — a„j + S — I. 



(30) 



V. CYCLOTRON RESONANCE OSCILLATOR 
STRENGTH 

The oscillator strength | < ^'■|'|pe''^|^' > p for electric 
dipole transitions between the LL is given by 



where 



T={Ta+Tb+ Ta' + Tb' 



7c = y 4>c,i^4>c,i'd^, 



(31) 



(32) 



with C = A,B,A' and B'. These matrix elements were 
calculated in the case of the Landau gauge, as 

Ta = didi, (^VlSi^i'+i + \/ITT^M'-i) / V2, (33) 



(a/ - 6){ai' - S) 



(34) 



Ta' = didi'fifi 



2V(/' + l)(; + l) 



Tp. 



{ai + S){ai> + 6) 
x(VTTT<5;,i'+r + VTT2V-i)/V2, (35) 

didi,fifi, (^VlSi^v+i + VITT(5m'-i) /V2. (36) 



where I and are Landau indices of the different eigen- 
states. The selection rule is SI = \l\ — \l'\ = ±1. Here 
(I, I') indicates a transition from a LL with index I to an- 
other with index V; the negative sign is used to represent 
a transition from a hole-like LL. The oscillator strengths 
for some dipole- allowed transitions between LL are shown 
in Fig. 6, for a biased bilayer (a) with Ui = —U2 = 50 
meV and i? = 10 T, as well as for an unbiased bilayer 
(b). The effect of the bias is most clear in the shift in 
energy of the (0,1) and (—1,-2) transitions, caused by 
the gap in the spectrum. On the other hand, one can 
also observe a decrease in the energies associated with 
transitions between the higher LL. 

Figure 7(a) shows the energy Ai? for some transitions 
between LL in a biased graphene bilayer as a function of 
the magnetic field, for Ui = —U2 = 50 meV, namely (0, 1) 
(black dashed lines), (1,2) (blue dotted line), (—1,-2) 
(black solid), and (—1,2) (red solid line). In contrast 
with the case of the unbiased bilayer as well as the single 
layer, AE is seen to be weakly dependent on the field 
for B < 5 T. For the (1, 2) and (-1, -2) transitions, the 
energy reaches negative values, due to the existence of 
LL crossings (see Fig 3). In addition, a difference in the 
energy of these transitions is found (see inset). The re- 
sults are obtained by subtracting the energies associated 
with the initial and final LLs. However, recent theo- 
retical work has indicated that many-body corrections 
may cause a substantial modification of the transition 
energies^^. The right panel Fig. 7(b) shows the oscilla- 
tor strengths for the same transitions as in the left panel, 
as a function of the magnetic field. These results show 
that, for small fields, e.g. for _B < 10 T, the oscillator 
strength associated with the (—1,2) transition is slightly 
larger than that of the (0, 1) transition, whereas for larger 
fields the latter becomes stronger. The curve for the 
(—1,2) transition intercepts the result for the (0, 1) one 
ai B 11.6 T, which corresponds to Ei — —Ui (see 
Fig. 3(b)). As in Fig. 7(a), the results show a significant 
asymmetry between the electron (1,2) and hole (-1,-2) 
intraband transitions. 
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This difference between tfie electron and liole intra- 
band transitions becomes larger as the gap parameter 
(Af7) is increased, as shown in Fig. 8, which displays 
the energies for the (1,2) (solid blue line) and (—1,-2) 
(dashed red line) transitions as a function of Af/ for 
B = 20 T, and (1,2) (solid green line) and (-1,-2) 
(dotted black line) for B = 10 T. At AU = 200 meV 
the energy difference reaches ~ 3 meV for B = 20 T, and 
« 0.8 meV for S = 10 T. For the B = 10 T case, the 
energy difference reaches zero and eventually becomes 
negative, due to the crossing of the LL. The asymmetry 
between the electron and hole intraband transitions is 
also particularly evident in their oscillator strength val- 
ues, as shown in Fig. 9, which presents the oscillator 
strengths for the transitions discussed in Fig. 7, as a 
function of the gap parameter, for B = 5 T. The results 
show that the presence of a gap causes a significant re- 
duction of the oscillator strength for the (0, 1) and (1, 2) 
transitions, whereas the result for (—1,2) shows a slight 
increase at small fields, reaching a maximum at AU ~ 25 
meV with a subsequent decrease. On the other hand, the 
curve for the (—1, —2) transition has a striking increase 
in comparison with the result for AU = 0. 

Figure 10(a) shows the energies for the (1,2) and 
(-1,-2) transitions (i.e. the cyclotron resonance tran- 
sitions) as a function of field, for several values of the 
gap parameter: AU = (black sohd line), AU = 25 
meV (red dashed lines), 50 meV (blue dotted hues), 100 
meV (green dot-dashed lines) and AU = 200 meV (ma- 
genta dot-dot-dashed lines). For AU = 0, the AE for 
the two transitions coincides, whereas in the other cases, 
the energy associated with the hole intraband transition 
is always higher than the corresponding electron tran- 
sition. Note that the AU = result only exhibit lin- 
ear _B-dependence, as predicted by the simple expression 
Eq. (27), for small _B-fields. With increasing AU the 
cyclotron energy decreases and for small i?- values AE 
becomes even negative, indicating a reversion of the or- 
der of the LL. These qualitative changes in AE should 
be observable in a cyclotron resonance experiment. 

The oscillator strengths for the transitions shown in 
the left panel are displayed in Fig. 10(b). It is seen 
that for the gapless system the oscillator strengths are 
independent of the magnetic field. On the other hand, as 
the gap increases, the results show a strong asymmetry 
that increases as AU becomes larger. In all cases, as 
B increases, the results approach the oscillator strengths 
for the AU ~ case. A similar behavior is also found 
for the other {I, I + 1) and (— /, —I — 1) transitions, albeit 
with smaller discrepancy in the energies. 

The asymmetry between the electron and hole intra- 
band transitions can be explained by taking into account 
the functional dependence of the amplitude factors of the 
spinor components on the energy. These terms, in the 
presence of a gap, introduce an asymmetry between the 
envelope functions of electrons and holes, thus modify- 
ing the transition probabilities. This effect is in turn a 
consequence of the breaking of the inversion symmetry 
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FIG. 6: (a) Oscillator strengths for electric dipole transitions 
in a biased bilayer of graphene {AU — 100 meV) for _B = 10 
T. (b) As in (a) for an unbiased ( AU = meV) bilayer. 



> 

0) 



160 




1 ' 




140 








120 




(-1,2)/ ; 


100 








80 


20 




(0,1) 


60 


> 






40 






(-1,-2) 


20 




B (T)" 


-^(1,2) ■ 







1 





10 15 
B(T) 



20 




5 10 15 20 
B(T) 



FIG. 7: (a) Transition energies in a biased {AU — 100 meV) 
graphene bilayer, as a function of the magnetic field B, for 
the dipole- allowed transitions (0, 1) (black dashed line), (1, 2) 
(blue dotted line), ( — 1,-2) (black solid line), and ( — 1,2) 
(red solid line), (b) Oscillator strengths vs field B for the 
transitions described in (a). 



of the biased bilayer. 



VI. SUMMARY AND CONCLUSIONS 

Using a four-band model we obtained analytical ex- 
pressions for the Landau levels and the eigenfunctions 
in a biased graphene bilayer in the presence of an ex- 
ternal perpendicular magnetic field B. In doing so 
we also extended previous results for the discrete spec- 
trum of an unbiased graphene layer based on a two-band 
model^^. Further, we obtained the selection rules for 
electric dipole transitions and evaluated the associated 
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FIG. 8: Transition energies as a function of AU for the (1,2) 
(blue solid line) and (-1,-2) (red dashed line) transitions for 
B = 20 T, and (1,2) (green solid line) and (-1,-2) (black dotted 
line) for B = 10 T. 
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FIG. 9: Oscillator strengths for dipole-allowed transitions in 
a graphene bilayer as a function of the interlayer potential 
difference AU for the (0,1) (black dashed line), (1,2) (blue 
dotted line), ( — 1, —2) (black solid) and (—1, 2) (red solid line) 
transitions, at B = 5 T. 



at the same time increasing the energy for transitions 
between an electron and a hole state, see Figs. 5 and 
6. Also noteworthy is the modification of the electronic 
band structure of the bilayer caused by the bias, which 
leads to the breaking of electron-hole symmetry, the ap- 
pearance of degeneracies in the LL spectrum, and quali- 
tative different dependency of the transition energy AE 
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FIG. 10: (a) Energies of the (1, 2) and ( — 1, —2) (upper curves 
of each doublet) transitions as a function of magnetic field for 
AU = (black solid line), AU = 25 meV (red dashed lines), 
50 meV (blue dotted lines), 100 meV (green dot-dashed lines) 
and AU = 200 meV (magenta dot-dot-dashed lines), (b) 
The oscillator strengths for the transitions shown in (a) as a 
function of magnetic field. The results below 0.92 correspond 
to the (1,2) transitions. 



on the magnetic field, which in turn can strongly mod- 
ify the dipole transition probabilities. These effects can 
be probed by cyclotron resonance measurements^, scan- 
ning tunneling spectroscopy^ as well as by far-infrared 
spectroscopy®. Future work may take into account the 
effects of spin as well as electron-electron interactions, 
which may become particularly important at high mag- 
netic fieldsZ'iiS'. 
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